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ABSTRACT 

We present the first detections by the NASA K2 Mission of oscillations in 
solar-type stars, using short-cadence data collected during K2 Campaign 1 (Cl). 
We understand the asteroseismic detection thresholds for Cl-like levels of photo¬ 
metric performance, and we can detect oscillations in subgiants having dom¬ 
inant oscillation frequencies around 1000 pHz. Changes to the operation of 
the hne-guidance sensors are expected to give signihcant improvements in the 
high-frequency performance from C3 onwards. A reduction in the excess high- 
frequency noise by a factor of two-and-a-half in amplitude would bring main- 
sequence stars with dominant oscillation frequencies as high as ~ 2500 /iHz into 
play as potential asteroseismic targets for K2. 

Subject headings: Astronomical instrumentation - K2 Mission 


1. Introduction 

Asteroseismology of solar-type stars has been one of the major successes of the NASA 
Kepler mission (Gilliland et ah 2010a). The nominal mission provided data of exquisite 
quality for unprecedented numbers of low-mass main-sequence stars and cool subgiants. 
Asteroseismic detections were made in more than 600 held stars (Chaplin et al. 2011a; 
2014), including a sample of Kepler planet hosts (Huber et al. 2013). These data have 
enabled a range of detailed asteroseismic studies (see Chaplin & Miglio 2013 and references 
therein), many of which are ongoing. 
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The nominal mission ended in 2013 May with the loss of a second of the spacecraft’s 
four onboard reaction wheels. This meant the spacecraft could no longer maintain three- 
axis stabilized pointing. However, thanks to the ingenuity of the mission teams, Kepler 
data collection has continued as a new ecliptic-plane mission, K2 (Howell et al. 2014). 
Targeting stars in the ecliptic minimizes the now unconstrained roll about the spacecraft 
boresight, thereby helping to compensate for the loss of full three-axis stability. The degraded 
photometric performance presents particular challenges for the detection of oscillations in 
solar-type stars. The oscillations are stochastically excited and intrinsically damped by 
near-surface convection. While this mechanism gives rise to a rich spectrum of potentially 
observable overtones, having periods of the order of minutes, it also limits the modes to tiny 
amplitudes, typically several parts-per-million in brightness. 

The opportunity to continue asteroseismic studies of solar-type stars with K2 would 
provide fresh data on stars in the solar neighborhood for application to both stellar and 
Galactic chemical evolution studies. The new helds have also led to the possibility of de¬ 
tecting oscillations of solar-type stars in open clusters and eclipsing binaries. This would 
provide independent data to test the accuracy of asteroseismic estimates of fundamental stel¬ 
lar properties. Other specihc targets of interest would potentially benefit from the provision 
of asteroseismic data, known exoplanet host stars being obvious examples. 

In this paper we report the detection of oscillations in several subgiants using K2 short- 
cadence (SC) data collected during Campaign 1 (Cl). We describe the target selection and 
data analysis, and also discuss the implications of our results for future K2 campaigns. 


2. Data 

2.1. Target selection and follow-up spectroscopic data 

Our selected target list started with the Hipparcos Catalog (van Leeuwen 2007). Use of 
these data allows us to make robust predictions for many bright, potential K2 targets in the 
ecliptic. Effective temperatures were estimated from the B — V color data in the catalog, 
using the calibration of Casagrande et al. (2010), and luminosities, L, were estimated from 
the parallaxes. These calculations used reddening estimates from Drimmel et al. (2003) 
(negligible for many of our targets). We adopted Mboi,© = 4.73mag (Torres 2010), and 
consistent bolometric corrections from the Flower (1996) polynomials presented in Torres 
(2010), which use the estimated Teg as input. We also applied a cut on parallax, selecting 
only those stars having fractional parallax uncertainties of 15 % or better. Stellar radii were 
then estimated from L and Tefj, and approximate masses were estimated from a simple power 
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law in L (which was sufficient for selecting targets). 

The estimated Hipparcos-based stellar properties were used as input to well-tested pro¬ 
cedures (Chaplin et ah 2011b) that enabled us to predict seismic parameters and relevant 
detectability metrics. We narrowed down the sample to 23 well-characterized bright (m^ ~ 7 
to 9) solar-type stars to be proposed for K2 observations. All targets were predicted to show 
solar-like oscillations on timescales of the order of minutes, necessitating SC observations. 

We also collected ground-based spectroscopic data on our selected Cl targets to help us 
check the Hipparcos-based predictions, and to better understand the hnal yield of asteroseis- 
mic detections. Observations were made using the TRES spectrograph (Fiiresz 2008) on the 
1.5-m Tillinghast telescope at the F. L. Whipple Observatory. Atmospheric parameters were 
derived using the Stellar Parameter Classihcation pipeline (SPC; see Buchhave et ah 2012). 
SPC was used to match observed spectra - taken at a resolution of 44000 - to sets of syn¬ 
thetic model spectra to derive estimates of Tes, log g, metallicity, and f sini. In what follows 
we assume that relative metal abundances [m/H] returned by SPC are equivalent to relative 
iron abundances, [Fe/H]. Table d] contains the derived spectroscopic parameters. There are 
four rapidly rotating stars in the sample, and some caution is advised regarding their esti¬ 
mated parameters. Overall, we found good agreement between the spectroscopic parameters 
and the Hipparcos-based values. Table [T] also includes the Hipparcos-based estimates of the 
luminosities. 

To understand the limits on K2 performance in Cl, we deliberately sampled the region 
of the HR diagram across which detections had been made in the nominal mission, as shown 
in the top panel of Fig. [H The symbols denote stars that provided hrm asteroseismic detec¬ 
tions (black), marginal detections (gray), no detections (open) or no detections with a high 
measured nsinz (red asterisks). Details are given below (notably in Section IXT]) . 


2.2. K2 Data and Lightcurve Preparation 

Each target was observed by K2 for just over 82 days in Cl—which lasted from 2014 
May 30 to 2014 August 21—with data collected in SC mode (Gilliland et ah 2010b). We used 
the K2P^ pipeline (Lund et ah 2015) to prepare SC lightcurves for asteroseismic analysis. 
In brief, the pipeline took the SC target pixel data as input. Masks for all targets in a given 
frame were dehned manually, and flux and position data were then extracted for our chosen 
targets of interest. Corrections were then applied to the lightcurves to mitigate the impact 
of changes of the target positions on the CCD. Finally, additional corrections were made 
using the hltering described by Handberg & Lund (2014). 
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Fig. 1.— Top panel: HR diagram of the SC targets. Targets with hrm asteroseismic de¬ 
tections are plotted with hlled black circles; marginal detections with hlled gray circles; no 
detections with open circles; and no detections with very high measured nsini (~ 38 to 
85kms“^) with red asterisks. Bottom panel: Symbols with errors show measured back¬ 
ground levels Rmax in the K2 spectra (in power spectral density units) at the predicted r'max, 
with rendering as per the top panel. Curves show threshold background power levels, Hthr, 
to have a ~ 90 % (black dashed) and ~ 50 % (black dotted) chance of making a detection. 
Uncertainties on the plotted thresholds are rendered in gray. Arrows show the result of re¬ 
ducing the excess high-frequency noise by a factor of two-and-a-half in amplitude (just over 
six in power). 
















































Table 1. Spectroscopic parameters and luminosities of targets 


EPIC 

HIP 

Teff (±77) 

log 9 (± 0 . 1 ) 

[Fe/H] (±0.1) Dsini (±0.5) Predicted 

log(L/L 0 ) 



(K) 

(dex) 

(dex) (kms“^) i^max (mHz) 

(dex) 


201162999 

56884 

5749 

4.44 

-0.13 

2.1 


-0.08 ± 0.07 

201164031 

56907 

5723 

4.47 

0.38 

0.9 


0.07 ± 0.05 

201182789^ 

57275 

6532 

4.15 

-0.04 

38.4 

1491 + 393 

0.69 ± 0.09 

201215315^ 

57456 

6523 

4.20 

-0.14 

41.3 


0.50 ± 0.04 

201343968 

55379 

6219 

4.09 

0.04 

9.0 

13491356 

0.61 ± 0.09 

201353392 

55288 

6110 

4.13 

-0.01 

8.1 

14651387 

0.43 ± 0.07 

201367296 

58093 

5695 

4.00 

0.19 

3.5 

113513™ 

0.42 ± 0.05 

201367904 

58191 

6125 

3.84 

-0.04 

10.3 

75411 ™ 

0.64 ± 0.11 

201421619 

55438 

5751 

4.39 

-0.38 

2.5 

27991740 

0.00 ± 0.12 

201436411 

56282 

6009 

4.15 

-0.26 

4.7 

1565l«? 

0.23 ± 0.06 

201592408 

56755 

5993 

4.31 

-0.24 

3.6 

226o1™7 

0.23 ± 0.11 

201601162 

54675 

5911 

4.53 

0.04 

5.6 

oQ-1 O+1007 
3812_.^gg 

0.07 ± 0.09 

201602813 

55022 

6156 

3.81 

-0.82 

11.9 

710HS 

0.33 ± 0.11 

201614568^ 

54857 

6940 

4.13 

-0.05 

85.8 

137813™ 

0.71 ± 0.03 

201620616 

58643 

5999 

4.37 

-0.16 

3.6 

2599l™0 

-0.02 ± 0.05 

201626704 

54541 

5505 

4.50 

0.09 

0.7 

36861™! 

-0.20 ± 0.03 

201698809 

55638 

5570 

4.71 

0.09 

0.4 

5916111 ™ 

-0.40 ± 0.05 

201729267 

55574 

6130 

3.76 

-0.43 

7.0 

622"*”^®^ 

*^^^-130 

0.80 ± 0.12 

201733406 

55467 

6155 

3.76 

-0.29 

5.3 

627"*”^®^ 

0.64 ± 0.09 

201756263^ 

57034 

6820 

3.81 

-0.08 

54.0 

674+17'^ 

0.56 ± 0.10 

201820830 

55778 

6417 

4.08 

0.00 

12.2 

1292 

0.62 ± 0.08 

201853942 

57136 

6053 

4.14 

-0.05 

4.4 

1 El 0 + 399 

iOi3_3j^Q 

0.33 ± 0.07 

201860743 

57676 

5852 

3.95 

-0.08 

4.7 

993+262 

^^'^ — 207 

0.56 ± 0.07 


Extra caution is advised regarding the classifications of these rapidly rotating stars. 








3. Results 


3.1. Asteroseismic Detections 


The lightcurves were distributed to several teams, who each attempted to detect sig¬ 
natures of solar-like oscillations in the power spectra of the data. A complementary range 


of well-tested analysis codes was used, which had been applied extensively to data from the 


nominal Kepler mission (e.g. Christensen-Dalsgaard et ah 2008; Huber et ah 2009; Mosser 
& Appourchaux 2009; Roxburgh 2009; Hekker et ah 2010; Kallinger et ah 2010; Mathur et 
ah 2010; Gilliland et ah 2011; Benomar et al. 2012; Campante 2012; see also comparison 
of methods in Verner et al. 2011). In cases where oscillations were detected, each team was 
asked to return estimates of the two most commonly used global or average asteroseismic 
parameters: Az/, the average frequency spacing between consecutive overtones of the same 
angular degree; and r'raax, the frequency at which the oscillations present their strongest 
observed amplitudes. 

We checked the asteroseismic detection yield using the spectroscopic data. The bottom 
panel of Fig. [T] provides a visual summary of these checks. The horizontal axis shows the 
predicted z/max for each target, made using the spectroscopic and logg as input. Esti¬ 
mates were calculated using the widely-used scaling relation (Brown et al. 1991, Kjeldsen & 
Bedding 1995): 



( 1 ) 


with the solar value z/max© = 3090/iHz (see Chaplin et al. 2014) providing the absolute 
calibration. These predicted z^max are given in Table [T] 

The vertical axis on the bottom panel of Fig. [1] relates to our ability to detect solar-like 
oscillations. In the frequency domain, where the analysis of the K2 data is conducted, peaks 
due to the oscillations are superimposed on a slowly varying, broad-band background. Our 
ability to make a detection depends on the prominence of the oscillation peaks above that 
underlying background. For solar-like oscillators, the background in the frequency range 
occupied by the most prominent modes has contributions from granulation, shot noise, and 
other instrumental noise. 

The symbols with error bars show measured background levels, Hmax, in the K2 spectra 
at the predicted r'max- Solar-like oscillators with lower frequencies of maximum variability 
show larger amplitude oscillations: the lower z /^ax , the larger is the maximum amplitude (i.e., 
brightness variation, in ppm) and hence the easier it is to make a detection for a given Rmax- 
Using the results from over 600 stars from the nominal mission we can predict oscillation 
amplitudes as a function of z/max- With the predicted amplitudes in hand, we may estimate 
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threshold background levels that would permit a signihcant detection of the oscillations. 

We calculated these thresholds using the detection recipe in Chaplin et al. (2011b). 
The lines in the bottom panel of Fig. [T] show the threshold levels, i^thr, below which the 
observed backgrounds i?max niust he to have a ~ 90 % (black dashed) and ~ 50 % (black 
dotted) chance of making a detection. The gray curves show uncertainties on the plotted 
thresholds. 

The results returned by the mode-detection teams indicated that we had four good 
asteroseismic detections. These stars are plotted in Fig. [T]using hlled black symbols; numbers 
are the hnal three digits of the associated K2 Ecliptic Plane Catalog (EPIC0) numbers (Huber 
et ah, in preparation). Two other targets showed marginal detections, and they are shown 
in gray. It is important to stress that these stars all lie in the part of the bottom panel where 
we would expect to make detections, i.e., where Hmax < -Bthr- Even though there are other 
targets that show lower or similar background levels in their K2 spectra, they he at higher 
predicted z/max where the intrinsic oscillation amplitudes, and hence the chances of making 
a detection, are lower. 

The observed z/max were in good agreement with the spectroscopic predictions, at the 
level of precision of the data. Note that the spectroscopic predictions of z/max have fractional 
uncertainties of ~ 25%, which are signihcantly larger than the typical uncertainties given 
by the asteroseismic measurements (which are in contrast at the few-percent level). 

Fig. [2] shows the K2 power spectra of the four stars with firm detections. Each spectrum 
has been smoothed with a 2 /iHz boxcar filter. Sets of vertical gray solid and dashed lines are 
separated by the estimated average Az/, and mark the spacing on which we would expect to 
see modes. The power envelope of the oscillation spectrum of EPIC 201820830 is somewhat 
flatter in frequency than the more classic Gaussian-like envelopes shown by the other three 
stars. This is a characteristic of hot F-type stars (e.g., Arentoft et al. 2008). The oscillations 
are strongly damped and this tends to wash out the visual appearance of the oscillation 
spectrum. 

The insets show the power spectrum of the power spectrum (PSPS) of each star, com¬ 
puted from the region around r'max- The significant peak in each PSPS lies at Az//2, and is 
the detected signature of the near-regular spacing of oscillation peaks in the frequency spec¬ 
trum. These detection signatures persisted for each star when the respective Cl lightcurves 
were divided into two equal halves in the time domain and analyzed separately. 

There are two stars in the the bottom panel of Fig. [1] that did not yield hrm or even 


^http://archive.stsci.edu/k2/epic.pdf 
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Fig. 2.— K2 power spectra of the four targets with hrm asteroseismic detections (smoothed 
with a 2 /xHz boxcar hlter). Insets show the PSPS of each star, computed for regions in the 
hrst spectra that lie about r'max- See text for further details. 
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marginal detections, but might be expected to do so. One of these stars (EPIC 201756263) 
is plotted in red. There are also three other stars shown in red on the diagram. These are 
all rapidly rotating solar-type stars with usini in the range ~ 38 to 85kms“^. They are 
therefore presumably young and may be very active, which is known to lead to signihcant 
suppression of the oscillation amplitudes (e.g., Garcia et ah 2010; Chaplin et ah 2011c; 
Campante et al. 2014) making it much harder to detect oscillations. The detection recipe 
does not yet make allowance for suppression of the oscillation amplitudes by activity. 

The other star that might be expected to show a detection (EPIC 201602813, plotted 
with an open symbol) is noteworthy because it is by far the most metal-poor star in the 
sample. This may lead to attenuation of the observed amplitudes, relative to the basic pre¬ 
dictions (e.g., Houdek et al. 1999; Samadi et al. 2010). Other explanations are that the z/max 
scaling underpredicts the z/max value for this star; or that the spectroscopic parameters are 
incorrect. We think that this latter explanation is unlikely. The target-selection (Hipparcos- 
based) prediction for z/^ax was in good agreement with the spectroscopic estimate. The 
target is also the primary component of a single-lined spectroscopic binary (Carney et al. 
2001), which will be spatially unresolved in the K2 pixel data. However, flux contamination 
from the secondary (a suspected white dwarf) is likely to be low. 

It is of course also possible that residual artifacts and other problems relating to the 
lightcurve extraction, preparation and hltering may have prevented a good detection being 
made. There are clearly persistent artifact peaks in some of the frequency spectra. Moreover, 
the high-frequency noise levels for all stars are not close to being shot-noise limited, unlike 
the nominal Kepler data. Analysis of the K2 SC spectra therefore demands some degree of 
careful, manual scrutiny to be sure that a claimed detection is not the result of a chance 
combination of noise peaks. 

Not including the active target EPIC 201756263, of the seven cases for which we would 
hope to detect oscillations, there are hrm detections in four and marginal detections in two. 
This is a good return given the challenges posed by the K2 photometry. 


3.2. Grid-Modeling Results 

We consolidated the Az/ and z/^ax estimates returned by the mode detection teams to 
give hnal asteroseismic parameters ready for asteroseismic modeling. The final parameters 
were those returned by an updated version of the OCTAVE pipeline (Chaplin et ah, in 
preparation; see also Hekker et al. 2010). Its estimates had the smallest average deviation 
from the median estimates in a global comparison made over all pipelines and all stars. 
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Uncertainties on each final parameter were given by adding (in qnadratnre): the formal 
parameter nncertainty given by the chosen pipeline; the standard deviation of the parameter 
estimates given by all other pipelines; and a small contribntion to acconnt for uncertainties 
in the solar reference values (which are employed in the grid-mo deling; see Chaplin et ah 
2014 for further discussion). The final parameters are listed in Table |2l 

A second set of teams were then asked to independently apply grid-based modeling 
to estimate fundamental properties of the four stars with firm detections. These teams, 
like their mode-detection counterparts, used codes that have been applied extensively to 
nominal-mission Kepler data (e.g., see: Stello et ah 2009; Basu et ah 2010; Kallinger et ah 
2010; Quirion et al. 2010; Gai et ah 2011; Bazot et ah 2012; Creevey et al. 2013; Hekker 
et al. 2013; Lundkvist et al. 2014; Miglio et al. 2013; Serenelli et al. 2013; Hekker & Ball 
2014; Rodrigues et al. 2014; Silva Aguirre et al. 2015). Further details may also be found 
in Chaplin et al. (2014) and Pinsonneault et al. (2014). 

We tested the impact on the estimated stellar properties of using different sets of inputs, 
i.e., a first set with {Au, z/max, [Fe/H]}, a second set with {Az/, Teg, [Fe/H]}, and further 
sets with the parallax-based luminosities also included. We found very good agreement 
between the properties given by the different sets, and by the different pipelines. Results 
from the first two sets were consistent with the independent luminosity estimates. Using 
the luminosities as an additional input constraint did not have a significant impact on the 
results. The parallax uncertainties—which range from ~ 5 % to 12 %—are too large to add 
anything useful to the seismic constraints. The fact that the second set of inputs provided 
results that were consistent with the estimated z/max lends further confidence to the claimed 
detections and suggests that, at least for these data, potential bias in the z^max estimates 
arising from spurious noise peaks is not a significant cause for concern. 

Table H] gives final values for the estimated properties, using the first set of inputs. 
The properties were calculated using the BeSPP pipeline (Serenelli et al. 2013), which uses 
individual model frequencies to calculate model predictions of Au for comparison with the 
observations. The uncertainties include a contribution from the scatter between pipelines 
(following the procedure outlined above for the input seismic parameters; see also the in- 
depth discussions in Chaplin et al. 2014). 

These grid-modeling results demonstrate that K2 has returned solid results that allow 
asteroseismic modeling to be performed on the targets. Two cases here - EPIC 201367296 
and EPIC 201367904 - will be amenable to more in-depth modelling studies since it will be 
possible to extract precise and robust individual frequencies of several overtones of each star. 
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4. Summary 

We analysed K2 short-cadence (SC) data for 23 solar-type stars observed in Cl. Of 
the seven targets where we would hope to detect oscillations, there are hrm asteroseismic 
detections in four cases, and marginal detections in a further two. This represents a good 
return, in spite of the challenges posed by the K2 photometry. In sum, we have a very good 
understanding of the asteroseismic yield. 

The results put us in a good position to hone target selections for future campaigns. 
Current performance levels mean we can detect oscillations in sub-giants, but not in main- 
sequence stars. Changes to the operation of the hne-guidance sensors are expected to give 
signihcant improvements in the high-frequency performance of K2 from C3 onwards. The 
high-frequency noise is currently a crucial limitation to making asteroseismic detections, in 
particular in main-sequence stars. With reference to the bottom panel of Fig. [1], we note 
that a reduction of excess high frequency noise by a factor of two-and-a-half in amplitude 
(just over six in power) would lead to consistent detections of oscillations in main sequence 
stars with z/max as high as ~ 2500 /iHz, as well as converting marginal detection cases to ones 
for which detailed modeling could be performed. 

The prospects are therefore very encouraging. Solar-type stars in the Pleiades and 
Hyades open clusters have already been observed by K2 in SC during C4. More stars will 
be observed in SC during C5 in the open clusters M44 and M67. There is also now clear 
potential to build up a statistical sample of solar-type held stars in the ecliptic that have 
good asteroseismic data, and to target specihc stars of interest for asteroseismic study such 
as bright eclipsing binaries and known exoplanet host stars. 


Funding for this Discovery mission is provided by NASA’s Science Mission Directorate. 
The authors wish to thank the entire Kepler team, without whom these results would not be 
possible. We also thank all funding councils and agencies that have supported the activities 
of KASC Working Group 1. 


Table 2. Asteroseismic parameters and estimated stellar properties 


EPIC 

HIP 

^max 

(mHz) 

Azz 

(mHz) 

M 

(M©) 

R 

(H©) 

P 

(gcm“^) 

log g 

(dex) 

201367296 

58093 

1176 ± 58 

65.7 ± 0.7 

1.14 ± 0.05 

1.71 ± 0.03 

0.323 ± 0.006 

4.032 ± 0.008 

201367904 

58191 

890 ± 46 

51.5 ± 1.0 

1.28 ± 0.05 

2.08 ± 0.04 

0.202 ± 0.006 

3.912 ± 0.009 

201820830 

55778 

1196 ± 72 

66.6 ± 0.8 

1.35 ± 0.06 

1.77 ± 0.03 

0.345 ± 0.006 

4.076 ± 0.008 

201860743 

57676 

1000 ± 46 

57.1 ± 1.3 

1.14 ± 0.05 

1.87 ± 0.04 

0.246 ± 0.007 

3.952 ± 0.010 
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